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Signals are a classical tool used in cellular automata constructions that proved to be useful for lan- 
guage recognition or firing-squad synchronisation. Particles and collisions formalize this idea one 
step further, describing regular nets of colliding signals. In the present paper, we investigate the use 
of particles and collisions for constructions involving an infinite number of interacting particles. We 
obtain a high-level construction for a new smallest intrinsically universal cellular automaton with 4 
states. 



Introduction 

Cellular automata were introduced by J. von Neumann flO'l in the 40s to study self-reproduction. They 
consist of a parallel computation model in discrete time based on an infinite grid of regular cells (endowed 
with a state chosen among a finite alphabet) interacting synchronously with each other. It is well known 
that even simple local interaction can lead to very complex global behavior. One open question is to find 
simple local function leading to a global rule, which is capable either of simulating any Turing machine 
or embedding any other cellular automata. Such automata are called respectively Turing or intrinsically 
universal. For a detailed survey of universality of cellular automata, see Q. 

For dimension 2 or above, the intrinsically universal cellular automaton constructed by E. R. Banks HI 
(2 states) is optimal. For dimension one, the result of M. Cook [2] (2 states) is optimal for the case of 
Turing universality whereas the previous best result on intrinsic universality is with 6 states [6|. 

All the previously mentioned constructions use structures known as particles and collisions. Those 
elements, which are commonly at the core of algorithmic on cellular automaton constructions, have 
been studied in-depth by J. Mazoyer and V. Terrier [5| . For constructions, a formal approach has been 
proposed in [8|. With this approach, it was possible to obtain [9| a new high-level proof of the result 
of M. Cook. In this paper, we use this formalism to construct quite simply an intrinsically universal 
cellular automaton with 4 states. The local transition function of the automaton is constructed to ensure 
the presence of particles and collisions with a predictable behavior. The technical part of the proof is to 
encode information with those elements to simulate any cellular automaton. Thus this proof is similar to 
complexity NP-completeness proofs. 

In section [T] we first recall definitions of cellular automata and different types of universality and 
then introduce the formalism used for particles and collisions. In section [2j we construct the rule of the 
cellular automaton and show how particles and collisions can be used. In section [3j we discuss in more 
details the encoding used and finish the proof. 
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1 Cellular automata, universality and self-organisation 

In this paper, we only consider cellular automata in dimension one, with the three nearest neighbours. 
Therefore, a cellular automaton is a pair {S,f) where 5 is a finite set of states and f : ^ S is the local 
transition function. The automaton acts on a configuration c S §^ at time t by synchronously applying 
the local transition function ;,c,-;,c,+i.f) to each cell c,- in c to obtain the value c^+i- The result of 

the global application of / to c is denoted by F{c). 

The automaton acts on c = (c,),^^ G (called configurations) by F(c), = /(c,_i , c,-, c;+i ). Moreover, 
a cellular automaton is one-way if the local function does not depend on its third argument. A space-time 
diagram D is, in our case, a bi-infinite sequence of configurations obtained in the dynamics of cellular 
automaton. More formally, it can be seen as an element of S^^ satisfying for all i,j G Z,D{i,j + 1) = 

f{D{i-lJ),D{i,j),D{i+^J))- 

Unlike many other known computation systems, cellular automata act on infinite configurations and 
do not posses any halting property. These two properties have led to several different definitions of 
universality. One position is to simulate Turing machine by using "ad-hoc" properties. A general scheme 
of different suggested properties was summarised by B. Durand and Zs. Roka |3|. Turing universality 
requires that, for each Turing machine, there exist a computable function that transforms each input word 
of the Turing machine into a ultimately periodic initial configuration and a pattern such that the initial 
configuration eventually containing the pattern if and only if the Turing machine eventually halts on the 
input word. 

The other approach is to simulate any other cellular automata evolution. Intuitively, a cellular au- 
tomaton simulates another cellular automaton if the space-time diagram of the latter can be "embedded" 
into the former. Formally, a cellular automaton 2t = {S,f) simulates a cellular automaton !B = {T,g) if 
there exists m,m' ,t,t' G N+, 5 G Z and a surjective encoding function e : S'" T'" such that, for any 
configuration c G T^, a'{F'{e-'^{c))) C e-\G'' (c)) where a : ^ is the shift defined for all c e 
and / G Nby a(c)(/) =c(/ + l). 

Definition 1. A cellular automaton is intrinsically universal if it can simulate any other cellular automa- 
ton. 

The notion of intrinsic universality is stronger than the notion of Turing universality in the sense that 
every intrinsic universal cellular automata is Turing universal but the converse is false. In the literature, 
two versions of intrinsic universality exist: one stronger form has the additional request that encoding be 
one-one. In fact most known constructions of universal cellular automata (except the one of M. Cook) 
achieve this stronger definition of universality. In this paper, we shall not enforce injectivity. Intuitively, 
this allows us to have some "junk" as long as it does not interfere with the simulation. Since we are 
interested mostly in the regularity of simulation, the two notions are equally interesting in our case. 
Furthermore, in a formal point of view, it is still not known whether they are the same. 

As many others constructions, we heavily rely on particles and collisions to encode and compute 
information. To define these objects, we use the approach developed in |8| seeing space-time diagrams 
as tilings of with local constraints. Therefore, we need to introduce some concepts from discrete 
geometry: a coloring is an application from a subset Sup{^) C 1? to S. Such coloring is finite if 
SupC^) is finite. Three natural operations on colorings are translation of a coloring by a vector u^I? 
defined by {u ■'^){z + u) = '^(z) ; disjoint union of two colorings and 'if' with Sup('^) n Sup('^') = 
defined by '^e'^'(z) = 'r(z)(resp. "^'(z)) for all z G Sup(<r')(resp. Sup(<^')); at last, restriction of a 
coloring to D C is denoted by "^/j. With those elements, let us give definitions of the three used 
structures of self-organisation. 
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(a) A background (b) A particle (c) A collision 

Figure 1 : Examples of structures present in self-organisation 

Structures are two-dimensional backgrounds (see Fig{T^) which are triplets *B = (^,m,v) where ^ 
is a finite coloring and u,v two non-coUinear vectors ensuring that 0,- j^j2{iu + jv) ■'rf is a. space-time 
diagram (this space-time diagram is often also referred as 55 when no confusion is possible). These 
structures mostly serve to "fill" empty space in constructions. Among backgrounds, one-dimensional 
particles travel (see Fig|I|)). Particles are quadruplets *p = {^,u,^i,^r) where is a finite coloring, 
u a vector, 03/ and !Br two backgrounds ensuring that ^ = ®kez^'^ ' separates the plane in two 4- 
connected domains L and R (L being the left-one according to u) so that !B|i © © iB'^ is a space-time 
diagram, where !B|i design the restriction of the coloring induced by *B on L. Particles are often used 
to convey part of information either alone or in groups consisting of parallel particles. Last but not 
least, a collision (see Fig. [T]:) is a pair ('^,L) where is a finite coloring, L is a finite sequence of n 
particles = (!B,-,'^-,m,-,!B-), satisfying the following conditions: first, consecutive particles on the list 
agree on their common background (for all / G Z„, 53- = 5B,+i), then particles and finite perturbation 
form a star ( = '^© 0,gz„ keN^'^i ' '^i ^^^^ the plane in n 4-connected zones and For all / G Z„, 
'^®®ken (^"i ■ '^i ® ^"(+1 ■ cuts the plane in two 4-connected zones. Let Pi be the one right of 
and at last, = © 0,- *B;|p. is a space-time diagram. 

To describe easily complex space-time diagrams, one idea is to symbolise particles as lines and 
collisions as points, giving birth to a planar map called catenation scheme as the one in upper-left corner 
of Fig. [3] Formally, a catenation scheme is a planar map whose vertices are labeled by collisions and 
edges by particles which are coherent with regards to collisions. Since catenation schemes are symbolic 
representations, it is not clear that there exist associated space-time diagrams. In fact, to go back from a 
catenation scheme to a "real" space-time diagram, one must give explicit relative positions of collisions 
as, for example, by giving the number of repetitions for each particle (edge) of the scheme. Such a 
set of integers is called affectation and is valid if the resulting object is a space-time diagram. The main 
result given for catenation schemes is that set of valid affectations can be computed from finite catenation 
schemes. 

Theorem 1 (f8l). Given a finite catenation scheme, the set of valid affectations is a computable semi- 
linear set. 

In this paper, the constructed automaton is based on particles and collisions and thus, heavily relies 
on the methodology used for catenation scheme. However, since the particles and collisions are explicitly 
constructed, they are very small and posses many good combinatorial properties. Therefore we do not 
need the whole power of catenation scheme and can often give simpler arguments for our specific case. 
The rest of the paper is devoted to construct the automaton and prove it is intrinsically universal. This is 
done in two steps: first, we give the automaton and show, using catenation scheme, that it can somehow 
"simulate" the local behavior of any automaton. Then, we show how to assemble those local simulations 
into a global one. 
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2 The automaton and elementary block 

In this section, we present the automaton and show how it can "somehow" simulate any local transition 
function of any one-way cellular automaton in an elementary block. Before going on with this block, one 
point to notice is that it is sufficient to simulate one-way cellular automata to be intrinsically universal. 
Therefore, our elementary block is constructed such that it takes three inputs (one transition function and 
two arguments) and outputs three values (one transition function and two "copies" of the result). 

To provide a good view of the automaton, all needed materials are depicted in Fig.[2]and Fig. [3] The 
first figure gives the local transition function and particles, whereas the second one gives the catenation 
scheme along with all interesting extracts of corresponding space-time diagram. The rest of this section 
is devoted to describe and explain the contents of these two figures. The local transition function is fully 
depicted on the top of Fig. [2] In fact, it is not a real transition rule since some cases (denoted by question 
marks) are not used and thus can take arbitrary values. 

Even if giving the rule could be sufficient, the next part is devoted to give intuitions behind local tran- 
sition function to ease understanding of the automaton. The first part is devoted to present background 
and particles used in the construction. Those structures are depicted in the corresponding part of Fig. [2] 
For each structure, we give a meaningful extract of space time, its name along with used local transition 
function cases. The formal definition can be trivially extracted from those diagrams. 

One requirement of our cellular automaton is that, contrary to other known constructions, it does not 
use a uniform background but a bi-colored check-board (B). With this new approach, we need two states 
(instead of one) for background but this allows us to have a greater range of particles since each remaining 
state can lead to two distinct particles according to its alignment within the background. Thus, with 
the two remaining states, we can construct four different particles (■,■,■,■). Furthermore, since 
background has two different phases, one can construct an additional particle (ff) by taking advantage of 
the gap between the two phases. Those are the main particles used in the construction. 

To encode information, the basic idea is to use groups of parallel particles that we call signals, 
depicted in the corresponding section of Fig. [2] In these groups, information can be encoded in two 
different ways: either by the number and type of the particles used or by relative position between those 
particles. Here, we use both approaches. As noted earlier, to simulate a on-way cellular automaton, 
we need to encode three main types of information: the left state (i.e., state of the left neighbour) is 
encoded in unary by the number of (regularly spaced) ■ particles (signal /). In a similar manner, the 
center state is encoded in unary by the number of (regularly spaced) ■ particles (signal m). The local 
transition function is encoded as an array of integers, with each integer encoded by spaces between a pair 
of particles ■ (signal R): the 7-th element of the array corresponds to the space between the j-th and 
the i + 1-th particle of the signal. To be exact, the space is counted by the sum of numbers of between 
two consecutive ■ and two (for the previously mentioned states). This way of counting may seem a little 
obscure but is chosen to give nice formulae in the end of this section. In the rest of the paper, names of 
symbols are also used when speaking of encoded values. 

During computation, other kind of signals appear: First, a mirror image of signal R (/?') which 
encodes the same information using ■ particles (going in the opposite direction). At last, an altered 
version of R' , which we denotes as R' , appears. In this altered version, some of the leftmost ■ particles 
are replaced by ■ . Due to choice of background and particles, our automaton has a special property 
which makes the proof a lot easier. Contrary to many other cellular automata, we have no synchronisation 
problems. 
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Local rule 




/(/,0,r) /(/.N,r) f{Lm,r) fUMr) 

Interrogation marks denotes cases that are not chosen and can be chosen arbitrary 



Background and Particles 




(0 (m) H^H (R) 

Collisions 



(a) ^ 'M\ \l7 ip9u 

(b) iM^ Jz zb 

(c) zk ^ ^ A Jz ^ 

(d) zb ^ \\- /?■-■■ ^ ^ bIz z^ 

(e) zk J- ^HB iHZ zSh aa^ Jz ^z 

(f) zkJ^AJz^AA^A 

(g) ^ 



Figure 2: Local rule and particles 
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Figure 3: Constructed block and collisions 
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Lemma 2. When two of the particles (or signals) described above collide, occurring collision is always 
the same. 

Proof. As studied by J. P Crutchfield et al. [T|, the number of ways two particles (or signals) can col- 
lide depends on the relative position between those two particles. Possible relative positions takes into 
account repetition vectors of those particles and constraints induced by backgrounds. In our case, repe- 
tition vectors are either (1, 1), (0,2) or (—1, 1) which suggests the possibility of two distinct collisions. 
However, the background forbids one of those possibility leaving only one possible case. □ 

With this very strong property, we can focus on symbolic behavior. Let us go and construct the 
dynamic. The scheme of elementary block is depicted in Fig. [3] (along with extracts of all induced 
collisions). The block is made the following way: left and right borders are delimited by particle. At 
the bottom, we have a left state signal /, a rule R and a center state signal m. The left state signal is 
going through the rule (collision /) and then collides with center state signal (collision c). This collision 
outputs an unused copy of the left signal to the right. This signal is erased when encountering the right 
border (collision g). Collision c also sends a signal encoding the the sum of left and center state to the 
left. This new signal is encountering (collision d) the mirror copy R' (created when R crossed the left 
border during collision b). During collision d, as many particle of R' are altere(|^ than encoded value 
(i.e., the sum of left and center signals). The signal embedding the sum is destroyed by the collision 
whereas the altered rule R proceed to the right. After crossing another rule R (collision e), the altered 
rule R collides with the right border and produces at the same time a new center state signal and a left 
state signal (located right of the border) in collision a. 

Additional cases needed for each collision in the local transition function are made explicit at the 
bottom of Fig. [3] This scheme gives us a symbolic block which somehow "computes" the rule. Now, let 
us prove that this symbolic behavior does really correspond to a valid space-time diagram and look at the 
details at the computed function. The following proposition deals with the first problem by ensuring that, 
under reasonable conditions, the scheme of symbolic block is a valid space-time diagram. Moreover, it 
gives some additional results on regularity of this block. 

Proposition 1. For any encoded rule R such that spaces between particles are even number greater that 
four, for any reasonable encoded value in left and middle state signal (i.e. both not null and such that 
their sum is less that the number of integers encoded in R) the scheme of the symbolic block correspond 
to a valid space-time diagram. Moreover, the size of the space time diagram does not depend on the 
encoded states and those blocks can tile the plane. 

Proof. First of all, the previous lemma ensures us that there is only one type of collision for each pair 
of particles; non null condition ensures that all particles exist. The proof that collisions are valid can 
be directly deduced from extracts given in Fig. [3] Due to the fact that constraints are local, periodicity 
considerations are sufficient to prove the validity of collisions occurring. In fact, it is sufficient to check 
the validity of one repetition of each periodic portion and the joint between different portion. The case 
of collisions c, e, f and g is easily get rid of since the perturbation inside the collision is periodic: adding 
one (or more) particles is the same as increasing the size of the periodic portion of the collision. Collision 
b requires space between particles ■ to be greater than four. For collision d, the constraint is just that 
there are at least one particles ■ left (i.e. sum of left and right values is less that the number of integers 
encoded in R). The last significant point is in collision a: since the vertical portion in the collision has 
periodicity (0,4), it requires that inter-space in R being even (this inter-space is of course the same as 

'Last alteration is an erasing rather than a change of the particle but this does not change the behavior 
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in R). If all these constraints are respected (which is the case for our proposition) then the resulting 
catenation scheme can be implemented as a valid space-time diagram. That is, there exists a space-time 
diagram where particles and collisions are positioned in the same way as in the catenation scheme. 

Now let us fix R, this implies that encoded values in states signals are bounded. Therefore, all signals 
are of bounded size and, up to increasing the size of the cell, they can be considered as objects with 
negligible size. With this, it is sufficient to take relative positions of signal as in the scheme of Fig.|3]to 
achieve same size blocks. 

The last point is to prove that such space-time diagrams, associated to elementary blocks, tile the 
plane. First remark is that symbolic blocks already tile the plane. Now if we look in details at Fig. [3] 
borders are left untouched by all collisions and rule signals are only shifted of a constant when encoun- 
tering a border. This implies that crossings of rule signals and border (collisions b) form a regular grid 
on the plane. At last, position of all other collisions only depends on a small number of neighbour points 
on this grid. It is directly visible for e and a. Positions of collisions /, g and c also depend on a' but its 
position is fixed by previous case. At last, collision d depends on c which was already treated. □ 

With this result, we have an elementary block able to do simple computation according to a rule R 
and which can tile the plane. Before combining those elementary blocks in the next section, we must first 
study exactly which function is computed by our block. For the same periodicity reason than previously, 
it is sufficient to look at what happens in the case of collision a in Fig. [3] 

Lemma 3. Let R(i) i<i<N be the array ofN integers encoded in R and x; ( resp. x,„) be the one encoded 
in I (resp. m); then the block leaves R unchanged and outputs m' and I' with encoded values respectively 
R{xi +x,n) +xi+Xm — {N + I) and R{xi +x,„) /2. 

This block computation is alike a cell in a cellular automaton except that it sends different values for 
left and right states. This difference prevents us to give a direct simulation and require quite additional 
work to get rid of this problem. A better way to overcome this problem would be to alter the automaton or 
use unused cases in the local transition function to ensure equality of outputs. However, for the moment, 
we do not manage to do such a thing. Therefore we present an alternative (and quite combinatorial) 
solution to use this elementary blocks in intrinsic simulation in the next section. 

3 Simulation of cellular automata 

To achieve intrinsic universality with our cellular automaton, the idea is to use the constructed elementary 
blocks and encode information not directly but only on portions of integers. First step is to remark that 
rather than working on only one block, we can work on chain of two consecutive blocks chained by left 
output (see Fig.|4]l. The new constructed element has three inputs (namely Xm) and three outputs 
xii, x,n' and x^' (not including rule). The key points is that we see input as integers written in binary and 
work on bits of these integers. The basic idea is to simulate a cell of a cellular automaton using x/ and 
Xni as inputs, x// and x„,' as outputs and maintaining the same constant value in all x,n and x,-/. 

Let us take a one-way cellular automaton {S,f). Since / does not depends on its third arguments, we 
can see it has an elements of — > S. Every element of S can be seen as a non null integer and written 
as a finite word of fixed size k on binary alphabet denoted by 5' = ^o^'i . . .Sk-i- All values encoded in 
signals can be seen as integers with 4 + 8^ + 3 bits (including leading zeros) as described in table [5] The 
first four one are called header, the last three one are called footer and each inside block of 8 (called 
digit) is used to encode one bit of the state. The size of the rule N is chosen as 2^+^^+^ — 1. Only some 
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Figure 4: Symbolic chaining of elementary blocks (encoding of state is omitted) 

parts of the signal for rule R encodes information. Other parts contain garbage which prevents a one-one 
encoding. 

The method used to encode and decode is depicted in Tab. [5] The encoding of / and m contains some 
junk but is chosen so that x/ +Xm contains the whole information about states encoded in / and m. With 
this, it is possible to chose the value of R{xi +Xm) according to these states. Half of the bits of /?(x; +x„,) 
are chosen to ensure the value of xp. The other half (denoted by symbol •) is used to ensure correct 
values in x^' = R{xi +x,„) +x/ +Xm — (A'^ + 1). In the second cell, one observe x;o +Xm which has a 
different header and thus a different value of any valid x/ +x,„. Therefore, one can chose /?(x/o +Xm) as 
N -l'^ + 1 (note that this value is even). With this choice, x,^,/ = N-f' + l + f + m-N- \ = x,n and 
X/' is on the correct form. With this encoding, we can simulate any rule of cellular automaton inside our 
cells. Leading to the following theorem: 
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Figure 5: Encoding of states inside values 
Theorem 4. The cellular automaton presented in Figgis intrinsically universal. 

Proof. With the previously presented encoding, chains of elementary blocks can encode any one-way 
cellular automaton local transition rule. Moreover, the encoding ensures conditions for proposition [1] 
Thus constructed elementary blocks can tile the plane. If we associate to each cell of the simulated cellu- 
lar automaton space-time diagram the corresponding (up to garbage consideration) chain of elementary 
blocks, we can simulate any other automaton. Last point is to note that even if the simulation is not 
horizontal, the problem can be easily overcome with standard trick. □ 



214 



A Particular Universal Cellular Automaton 



Conclusion 

Our construction heavily relies on signals and manages to give a construction with two separate levels: 
local rule is constructed to ensure a set of particles and collisions. Then, the simulation is made encoding 

the computation with this set. This method allows us to have a clear and understandable construction 
which separates the local and global aspect and gives us a new smallest intrinsically universal cellular 
automaton. 

References 

[1] E. R. Banks (1970): Universality in cellular automata. In: Symposium on Switching and Automata Ttieory 

(Santa Monica, California, 1970). ffiEE, pp. 194-215. 
[2] M. Cook (2004): Universality in Elementary Cellular Automata. Complex Systems 15, pp. 1^0. 
[3] B. Durand & Zs. Roka (1999): The game of life: universality revisited. In: M Delorme & J Mazoyer, editors: 

Cellular automata: a parallel model. Kluwer, Dordrecht, pp. 5 1-74. 
[4] W. Hordijk, C. R. Shalizi & J. R Crutchfield (2001): Upper bound on the products of particle interactions in 

cellular automata. Phys. D 154(3-4), pp. 240-258. 
[5] J. Mazoyer & V. Terrier (1999): Signals in one-dimensional cellular automata. Theoretical Computer Science 

217(1), pp. 53-80. Cellular automata (Milan, 1996). 
[6] N. OUinger (2002): The Quest for Small Universal Cellular Automata. In: ICALP '02: Proceedings of the 

29th International Colloquium on Automata, Languages and Programming. Springer- Verlag, London, UK, 

pp. 318-329. 

[7] N. OlUnger (2008): Universalities in Cellular automata; a (short) survey. In: B. Durand, editor: Proceedings 
of the First Symposium on Cellular Automata "Joumees Automates Cellulaires". pp. 102-118. 

[8] N. OUinger & G. Richard (2008): Collisions and their Catenations: Ultimately Periodic Tilings of the Plane. 
In: IFIP-TCS, 273. Springer Boston, pp. 229-240. 

[9] G. Richard (2008): Rule 110: universality and catenations. In: B. Durand, editor: Proceedings of the First 
Symposium on Cellular Automata "Journees Automates Cellulaires". pp. 141-160. 
[10] J. Von Neumann (1966): Theory of Self-Reproducing Automata. University of Illinois Press, Urbana, 111. 



